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Abstract 

White noise-driven nonlinear stochastic partial differential equations (SPDEs) of parabolic 
type are frequently used to model physical systems in space dimensions d — 1,2,3. Whereas 
existence and uniqueness of weak solutions to these equations are well established in one 
dimension, the situation is different for d > 2. Despite their popularity in the applied 
sciences, higher dimensional versions of these SPDE models are generally assumed to be ill- 
posed by the mathematics community. We study this discrepancy on the specific example 
of the two dimensional Allen-Cahn equation driven by additive white noise. Since it is 
unclear how to define the notion of a weak solution to this equation, we regularize the noise 
and introduce a family of approximations. Based on heuristic arguments and numerical 
experiments, we conjecture that these approximations exhibit divergent behavior in the 
continuum limit. The results strongly suggest that shrinking the mesh size in simulations of 
the two-dimensional white noise-driven Allen-Cahn equation does not lead to the recovery 
of a physically meaningful limit. 

Keywords: stochastic PDEs; well-posedness; Allen-Cahn equation; white noise; numerical 
analysis 



1. Introduction 

Stochastic equations of the type 

dtu = Au + f(u) + £, (1) 

where A is a linear elliptic differential operator, / is a nonlinear function and £ is space-time 
white noise, play a central role in the modeling of a whole variety of phenomena in the phys- 
ical sciences. Prominent examples are the Swift- Hohenberg equation in the study of thermal 
convection [44j or the Kardar-Parisi-Zhang (KPZ) and Lai-Das Sarma- Villain equations in 
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surface growth [27], [3JJ] . Another important domain of application for nonlinear parabolic 
SPDEs driven by additive white noise is the theory of the dynamics of critical phenomena 



18l |4|, [24J . Hohenberg and Halperin [24| introduced a nowadays widely used classification 
of the various models of dynamic critical phenomena. Their classification includes among 
others: Model A for non-conserved quantities, e.g. the stochastic time-dependent Ginzburg- 
Landau and Allen-Cahn equations; Model B for conserved quantities, e.g. the stochastic 
Cahn-Hilliard equation; Model C which couples conserved and non-conserved fields, e.g. 
phase-field models of eutectic growth. 



Over the past decade, the applied science community has paid a lot of attention to the 
two- and three-dimensional versions of these white noise-driven models. In addition to some 
analytic work, see e.g. [48j for the KPZ equation, emphasis has been put on numerical inves- 
tigations: standard finite difference and pseudospectral methods for spatial discretization 
have been combined with Euler-Maryama or stochastic Runge-Kutta schemes to find numer- 
ical approximations to these SPDEs. Fo r sp ecific examples of such work we refer to 
for the KPZ equation, to @, 0, [H, [35 , 
for Model C. 



for Models A and B, and to [H, [13, E3 



33. 
12, 



17] 
11] 



But while these models are extensively studied by applied scientists in space dimensions 
d — 1,2,3, in the case of nonlinear parabolic SPDEs with additive space-time white noise 
the mathematical literature focuses almost exclusively on the case d = 1. In fact, in one space 
dimension, the theory of nonlinear parabolic SPDEs of generic type (H} is well-established 
in the literature, see e.g. @, 45|. In addition, convergence properties of standard numeri- 
cal approximation techniques for such equations in K 1 have been thoroughly studied, see 



22l [2lJ, [32, |4l| for finite difference and [46[ for finite element methods, respectively. On the 
other hand, little analytic work has been done on the higher-dimensional cases of these white 
noise-driven models. With the exception of a few specific cases such as the stochastic Cahn- 
Hilliard equation (whose well-posedness is established in [(|), most mathematical studies 
claim upfront that additive white noise leads to ill-posed equations in d > 2. Typically, 
authors then follow one of the following regularization procedures: they either restrict the 
analysis to the case of colored noise with a finite spatial correlation length 2t| T3l [28l l26l. [30jj . 
or they render the equation well-posed by means of the so-called stochastic quantization pro- 
cedure involving Wick products @, Q • 



Assembling the above observations, we end up with the following discrepancy: whereas 
higher-dimensional versions of these nonlinear parabolic SPDE models with additive white 
noise are commonly used by applied scientists, they are assumed to be ill-posed by the 
mathematical community. The goal of this study is to gain a better understanding of what 
goes wrong in a class of Model A equations, and we do so by focusing on the specific example 
of the stochastic Allen-Cahn equation on a bounded domain in R 2 , 

d t u = Au + u - u 3 + f , (2) 

where £ is space-time white noise. As discussed in Section 13.11 below, it is unclear how 
to define a weak solution to this equation. Therefore we take a different approach: we 
consider a sequence of regularized versions of ([2]) with corresponding solutions un, and 
study the limit lim;v->oo u^. In particular, the regularization procedure is such that the 
sequence {wjv(t)}^ =1 is £ 2 -valued for all t > 0. As illustrated in Fig. [TA-D, the sequence 
of regularizations (shown at time t = 1) does not seem to converge to a meaningful limit as 
N — >• 00: the underlying deterministic evolution of the field (Fig. [TJC) gets washed out as 
N increases and the field becomes highly oscillatory (Fig. [Tp-D). Furthermore, the integral 
of ujv(i) over an arbitrary subset of the spatial domain decays with increasing N (Fig. [If). 
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Figure 1: 2D stochastic Allen-Cahn equation. A-D: Equation 1 160 I I is solved over time interval [0, 1] 
using scheme ifST) for increasing number of grid points: N = 8 (A), N = 32 (B), N = 128 (C), N = 512 (D). 
The final fields u^(x, y) are shown from side and top angles; values outside of [—5, +5] are set to —5 and +5 
respectively. Initial condition uo(x, y) = sin(2x); parameters: a = 6.4 ■ 10 -3 , g = 0.5, a = 2 -3 tt, M = 10 3 . 
E: The deterministic Allen-Cahn equation is solved using scheme 1611 . Initial condition and parameters 
as above, except a = 0, TV = 256. F: [— tt,+tt] is divided into four subintervals Iy. = (k — 1), ^k] for 
k = — 1, . . . , 2. Simulations A-D are repeated for 120 realizations of the noise and means of the piecewise 
constant functions «jv(z £ Ik) '■= fx Jq" u m (x , y) dy dx are plotted. For each 1^, the largest errorbar is 
shown. 

This motivates the following conjecture. 

Conjecture 1.1. The regularizations un of the two-dimensional stochastic Allen-Cahn 
equation converge in probability to the zero-distribution, i.e. Vt > and all smooth 
test functions </>: 

lira (un (t) , <f>) — > in probability, (3) 
where (•, •) denotes the L 2 - duality pairing. 

A simple model of this form of convergence consists of the sequence of functions {sin(A r x)}^ =] 
with x € [0, 2ir]. All elements of this sequence are smooth functions, but they become highly 
oscillatory for large N, similarly to Fig. [Tp. It is easy to show that for all periodic and 
smooth test functions <j), 

lim (sin(AT-),0) = 0, (4) 

N— too 
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which means that sin(iV:r) weakly converges to zero. 

The main goal of this study is to provide sound numerical evidence for Conjecture 11.11 and to 
give a heuristic explanation as to why the limit lim^^co un is not meaningful. After a brief 
review of the deterministic Allen-Cahn equation and the properties of its stochastic version 
for d = 1 in Section [5J we proceed to two simplified versions of the equation in general 
dimension d > 1. In Section [3] we revisit the well-studied stochastic heat equation and in 
Section2]we consider the decoupled stochastic Allen-Cahn equation - a field of independent, 
noise-driven particles in double-well potentials. Studying these simplifications provides us 
with a better understanding of the respective roles played by the nonlincarity and the 
diffusion in the full stochastic Allen-Cahn equation in two space dimensions. The latter 
is the subject of Section [5j where we provide heuristic arguments and numerical evidence 
for Conjecture 11.11 In Section [6] we present the numerical solution strategies employed in 
Sections [3][S] Finally, implications of our work for published numerical studies as well as a 
presentation of future directions are the subject of Section [7] 



2. Preliminaries 

The Allen-Cahn equation on a domain il C K d , 

d t u = Au + i (u-u 3 ) , (5) 

was first introduced by Allen and Cahn to describe a non-conserved order field during 
anti-phase domain coarsening [l[. The equation describes the gradient flow of the energy 
functional 

E[u] = [ ~ \\7u\ 2 - \v(u)dx (6) 

with the particular choice of a double- well potential V(u) = ju 4 — \u? . The Allen-Cahn 
equation has been extensively studied in the literature and exhibits a variety of interesting 
phenomena such as interface motion by mean curvature flow in the limit as e — > 0, see [14j j . 
We are not concerned with such limits and shall hereafter set e = 1. 



To take into account thermal effects in the gradient flow of ((SJl, a natural approach is to 
complement the flow with a random forcing term in form of additive space-time white noise 

d t u = Au + it — u 3 + cr£, (7) 

where a > is a constant and £ is a space-time Gaussian random process with mean zero 
and formal correlation function E t)£(x' , t')] — 8(x — x') 8(t — t') 42]. Note that in 
the framework of Hohcnbcrg and Halperin's classification [24|, the stochastic Allen-Cahn 
equation ([7]) belongs to Model A: systems without conservation laws. In the current study, 
we consider ([7]) on the torus T d , i.e. the d-dimensional hypercube [— 7r, +ir] d subject to 
periodic boundary conditions. This yields the following initial value problem on T d : 

( d t u = Au + (it - u 3 ) + of, t>0 

\ u = u , t = 0. {) 

For the sake of simplicity we assume that uq is a deterministic function. Later on, the 
following concise definition of space-time white noise for the torus-setting will be useful. 
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Definition 1. Let T > and consider the measure space ([0,T] x T d ,B, A), where B is the 
Borel a-algebra of [0,T] x T d and A i/ie Lebesgue measure. Let (fi,.F, P) fee a probability 
space. A space-time white noise on [0,T] x T d is a mapping (:6^ L 2 (Vt) such that 

(i) for all B € B, is centred Gaussian with 

E [£(B)] 2 = \(B). 

(ii) if Bx n ... n i? n = 0, £/ien £/ie are independent and 

n 

e(us =1 B i ) = E^i). 

Recall now that white noise is, formally speaking, the time derivative of the infinite di- 
mensional cylindrical Wiener process W [37], Sect. 7.1.2]. We can thus rewrite problem (|8]) 

as 

J du = [Au + u- u 3 ] dt + adW, t>0 , . 

\ u = u , t = 0. ( > 

This notation will be convenient for our analysis because W admits the spectral decompo- 
sition @ 

W(t) = Pk(t)e k , (10) 

where {/3fc} fcgZ d are i.i.d. Brownian motions and {ek} kGli d is an orthonormal basis of L 2 (T d ) 
with respect to the inner product (h, g) = J"r hgdx, 

e k (x) = (27r)- d/2 e* kx . (11) 

Before attempting to solve ((8]), we have to define the notion of a solution. In fact, white 
noise is too rough to make sense of the equation pointwise, and we proceed formally by 
integrating ([§]) against a smooth test function to obtain the following weak formulation [9( . 

Definition 2. Let % be a Hilbert space. An %-valued process u{t), t G [0, T], is called a 
weak solution to problem |5P if (i) J \\u(t)\\dt < +oo for almost all trajectories, and (ii) 
it satisfies the weak formulation 

(u(t), 0) = («o, 4>) + [ [<«(«), A0> + (u(s) - u 3 (s), 0)] ds + (W(t), </>) , (12) 
Jo 

almost surely, for all t € [0, T] and for all (f> G C°°(T d ). Here, (•, •) denotes the inner product 
on H. 

For the case d = 1 , existence and uniqueness of the weak solution (fT2j) have been established 
for U = L 2 (T) [Hi ( in fact, the solution is almost surely continuous). Therefore, the 
initial value problem ((5J) provides a mathematically sound model in one space dimension. 
However, in higher dimensions (d > 2) the situation is quite different: the weak solution to 
the linearized version is only a distribution-valued process, and hence it is unclear whether 
there is a Hilbert space % such that a unique weak solution to the nonlinear equation exists. 
To gain a better understanding of these issues, it is instructive to have a closer look at the 
stochastic heat equation, a linear version of the full Allen-Cahn model. 
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3. Simplified version I: stochastic heat equation 

3.1. Analytic Considerations 

In this section we focus on the well-studied stochastic heat equation on T d [iBj ]. 
subject to homogenous initial conditions, 

du = [Au -u]dt + adW, t > . . 

u = 0, t = 0. ( ' 

Using the spectral decomposition of the noise (jTUJ) , the projection of (TlU)) onto the Fourier 
modes (flT]) yields the following set of stochastic differential equations for k £ 1 d 

du k = -fikUkdt + ad/3 k , t>0 . . 

u fc = 0, i = 0. 1 > 



where fik = 1 + |fc| 2 and |fc| 2 = 5Z,- =1 In other words^the solution to (ITBl is represented 
by a collection of i.i.d. Ornstein-Uhlenbeck processes ([16], pl06) whose solutions are given 

by 



fi fc (f) = a / e-^( 4 - s )d/3 fc ( S ). (15) 
Jo 

These are mean zero Gaussian processes with covariance (s > 0) 

9 

Eu k (t)u k {t + s) = ^—e~^ s [1 - e" 2 ^] . (16) 

With this we are now able to calculate the expected value of the L 2 (T d )-norm (denoted by 
|| • ||o) of the solution to (fT3j) . 

E\\u(t)\\l = E £ \u k (t)\ 2 = £ ^ t 1 - e ~ W ] ■ ( 17 ) 

The convergence of the series depends on the summability of /i^ 1 = (1 + |fc| 2 ) -1 and it is 
easy to see that E||it(£)|| 2 < +oo for d = 1, whereas E||m(£) || 2 = +oo for d > 2. In other 
words, the solution of the one dimensional stochastic heat equation ()13[) takes values in 
L 2 (T X ) almost surely - a result that does not hold true in higher dimensions. 

Let us now have a closer look at the case d>2, and more precisely at the rate of divergence 
of the sum in ([P?). For fixed t > 0, 

i /•w d-l 

nnml- £ ?-^W rr^ dr ' (18) 

1^1 + 1*1 -/o i + T- 

where we use the following definition: /jv ~ ffjv if there exist two constants c, C > such 
that c < /at /<7jv < C for iV sufficiently large. (fTS]) implies that the divergence is logarithmic 
for d = 2 and polynomial for d > 3. But even though it is not L 2 (T d )-valued, u(t) might 
be well-defined in a larger space. In this sense, a natural enlargement of L 2 (T d ) is given by 
the Sobolev spaces of negative order, H s (T d ) for s < (e.g. [10(, p96). 

Definition 3 (Sobolev Spaces). Let set. Then the Sobolev space H s (T d ) is defined by 
means of a weighted norm as 



H°(T d ):=lf:T d ^R : ||/|| 2 = £ (l + |fc| 2 ) f k < +oo \ , (19) 
where fk are the generalized Fourier coefficients. In particular, L 2 (T d ) = H°(T d ). 
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Revisiting the calculation (fT5|) in H s (T d ) yields 



fN d-1 

E ll<*)lla~y o (1+r2)1 _> <+oo, Vs<l-d/2, (20) 

i.e. the solution to the <i-dimensional version of (fT5|) takes values in H s (T d ) almost surely, 
for all s < 1 — d/2. We summarize these estimates in the following result [45| . 

Theorem 3.1. Let d > 1. for a/Z s < 1 — d/2, the solution to the d-dimensional stochastic 
heat equation takes values in H s (T d ) almost surely. In particular, d — 2 is the borderline 
case: the L 2 -divergence in H8]) is logarithmic, and for all t > we find that u{t) E H s (T d ) 
almost surely, for all s < 0. 

Theorem 13.11 illustrates the smoothing properties of the heat kernel: whereas space-time 
white noise only takes values in H s (T d ) for s < —d/2, the action of the heat semigroup 
improves the regularity such that u(t) E H s (T d ) for s < —d/2 + 1. In view of the nonlinear 
stochastic Allen-Cahn equation ([8]) in dimensions d > 2, the relevance of Proposition 13.11 
is the following: for negative s, H s (T d ) is a space of distributions, and it is in general not 
possible to multiply two distributions in a meaningful way [40] . In consequence, taking the 
cube of the linear solution does not make sense and we anticipate the nonlinear equation to 
be ill-posed. This observation provides a heuristic argument as to why white noise-driven 
nonlinear parabolic SPDEs in higher dimensions are generally suspected to be ill-posed. 
However, it is not a rigorous proof and the question remains as to whether it is possible to 
find a space in which there exists a weak solution to the stochastic Allen-Cahn equation in 
higher dimensions. 

3.2. Simulations 

To illustrate the results of the previous section, and to facilitate comparisons with 
subsequent results, we present now a numerical experiment on the two-dimensional stochas- 
tic heat equation (JT3J). We discretize the periodic domain T 2 using A 2 grid points {xj = 
t5 j : j — j'2), ji — —N/2, . . . , N/2 — 1} and denote by ujy(xj,t) the numerical approx- 
imation at grid-point Xj and time t. Since we are mostly interested in the spatial regularity 
of the solution, we fix T > 0, integrate the equation numerically over [0,T], and denote 
Ujv(xj) = un(xj,T). To estimate the regularity of the approximation ujv, we first need to 
define the finite-dimensional analogue of the Sobolev norm (fH))) . To do this, we introduce 
the discrete Fourier transform of as 

1 N/2-1 N/2-1 

UN(k) = — u N(x 3 )exp(—k-j\, (21) 

h=-N/2 h=-N/2 ^ ' 

and the discrete inverse Fourier transform as 

J N/2 — 1 N/2-1 ^ 

u n(x 3 ) = — il N(k)aq>(- — k-j\, (22) 

k 1 = -N/2 k 2 = -N/2 ^ ' 

where {k — (ki, k 2 ) : fc, = —N/2, . . . , N/2 — 1} are the wave vectors. We can now define 
the iJ s (T d )-norm for discrete functions defined on the A 2 grid points. 

Definition 4 (Discrete Sobolev Norm). Let N > 2 be an even integer. Let be a discrete 
function defined on the N 2 grid points 0/T 2 . The discrete H s (T 2 )-norm of is defined 
as 

N/2-1 N/2-1 

\\\u N \t~p- 1 E (l + |fc| 2 ) S |^(fc)| 2 , (23) 

k 1 = -N/2 k 2 = -N/2 
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where p 1 ' 2 := 2n/N is the grid spacing. 

Remark 1. To see that p~ x is the correct scaling in \23\) , consider the case s — 0. Using 
Parseval's equality we find 

N/2-1 N/2-1 N/2-1 N/2-1 

lllwjvlllo = P~* Y Y \u N (k)\ 2 = p^ 1 Y Y \ u n(x)\ 2 

k 1 = -N/2 k 2 = -N/2 x 1= -N/2 x 2 = -N/2 

^=^> / \u\ 2 dx=\\u\\l (24) 

J[0,2tv] 2 

The following notion of the radial energy density will be useful for the graphical illustration 
of the numerics. 

Definition 5 (Radial Energy Density in Fourier Space) . Let N > 2 be an even integer and 
let un be a discrete function defined on the N 2 grid points of T 2 . The radial energy density 
in Fourier space, Em , is defined as 

E n{ k ) : = Tj-r Y l"Jv( fc )| 2 5 k = 1,..., N/2-1, (25) 

where A K = |a; G K 2 : (k — l) 2 < \x\ 2 < k 2 | is the n-th annulus in R 2 and \A K \ is the 
cardinality of 7? PI A K . 

Using the energy density (|25)l . we can replace the 7I s (T 2 )-norm (l23l) with an equivalent 
discrete norm defined as 

N/2-1 

lll«Jv|||»= Y kE n (k)(1 + k 2 ) s . (26) 

K = l 

From this we conclude that in the continuum limit as N — > oo, the regularity of is 
determined by the decay of the radial energy density Em- In Figure [5J the density Em for the 
corresponding numerical solution um of the stochastic heat equation is plotted for increasing 
values of N: we observe convergence within the errorbars, and in particular the Em(k) 
decay like 1/k 2 for k sufficiently large. By (l26l) . this suggests that limjv-j-oo IE||itjv||o = +°°; 
whereas limjv-^oo um € H S (T 2 ) almost surely, Vs < 0. Note that these limits are in perfect 
agreement with Theorem 13.11 A detailed discussion of the employed numerical scheme is 
presented in Section HTT1 

4. Simplified version II: decoupled stochastic Allen-Cahn equation 

4-1- Analytic Considerations 

In this section we study a second simplification of the stochastic Allen-Cahn equa- 
tion in W l . This time, we keep the nonlinearity but neglect the diffusion term, i.e. we 
consider 

du=[u- u 3 ] dt + adW, (27) 

on the domain T d and subject to homogenous initial conditions. Even though this equation 
does not provide a meaningful physical model, it is instructive in view of the discussion 
of the full stochastic Allen-Cahn equation in Section [5] In fact, we will see that approxi- 
mations to the stochastic Allen-Cahn equation and its decoupled version (|2"7) share similar 



8 



LU 



10 

10"" 
10~ 
10" 
1Cf 

io J 

10" 

io- 

10"' 




10 



10 



10 



10 



Figure 2: 2D stochastic heat equation. Equation H55II is solved using scheme H57I) , Parameters: g = 1, 
a = 0.5, a = tt/50, T = 1, M = 2 ■ 10 3 . For each N = 2" (n = 5, . . . , 11), the sample mean of E n (k) 
over 40 simulations is plotted for k = 1, . . . ,N/2. Largest error bar is shown for each N, and we observe 
convergence within errorbars. Dotted line is of slope 1/k 2 . 



characteristics. However, the latter is more tractable due to the spatial decoupling: at each 
point in space there is a particle, driven by a Brownian motion and confined by a double- 
well potential. Note that the particles are i.i.d. because the collection of driving Brownian 
motions arises from space-time white noise. It is not clear if there exists a solution space 
which allows for the definition of a weak solution to (|2"T)) . We circumvent this difficulty 
by first discretizing the equation in direct space and then passing to the continuum limit. 
We choose an even integer N, introduce an equidistant mesh of size p -1 / 2 := Ax = 2-k/N 
along each dimension of T d , and consider the finite-dimensional approximation un, which 
is defined on the N d grid points {xj}. At each grid point, the evolution equation is 



du N (xj) = -V {u N (xj)) dt + op d/A dfi 



(28) 



where V{x) :— double-well potential, {Pj} is a collection of i.i.d. standard 

Brownian motions. The scaling p d ^ 4 is due to the spatial discretization of white noise. 
The mesh partitions the domain into N d hyper-cubes of volume p~ d / 2 , and from Definition 
[1] it follows that the average noise on each hyper-cube (at fixed time i) is distributed as 
p d / 2 Af(0, p~ d l 2 ) ~ p d / 4 A/"(0, 1). As the number of grid points increases, the deterministic 
term in (|28p remains unaltered while the noise intensity increases. Therefore, the variance of 
the solution ujv(i) is unbounded in the limit as N — > oo. Unless the potential V is quadratic, 
it is not possible to write down the dynamic solution to (|28[) in closed form; instead we focus 
on the stationary solution (as t ~ > oo) whose probability distribution function is given by 



Pn(v) = ^ ex P 



a 2 p d/2 



(29) 



with J\f the normalization constant [16| . We prove now the following result 
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Theorem 4.1. Let d > 1 and V(x) = ^ir 4 — \x 2 . Then there exists a constant K > such 
that the stationary solution un of the regularized problem i28\) admits the following limits 

if s < -d/4 

lim E |||uiv|||a = i K < +oo if s = -d/4 (30) 
+oo if s > -d/4, 



where |||-||| s is the discrete Sobolev norm of Definition^ generalized to d dimensions A35\) . 

Proof. The main part of the proof is established in Lemma |H] in the Appendix. We first 
generalize the discrete Fourier transform (|2"TT) and its inverse (|2"21 to d dimensions: 

N/2-1 . . . 

^N{k) = — rF r V u N (xj)exp — k-j , (31) 



and the inverse transform 



N d/2 »JV V ~J/~V I Ar 

3i=-N/2 



N/2-1 . „ . 

u N(xj) = u N{k)exp(-—k-jj. (32) 

k i = -N/2 ^ ' 

Using (f3~T|) together with point (ii) of Lemma [9j we find that the stationary solution un of 
(|2"8)l satisfies the following properties: 

EuAr(fc)=0, Vk = -N/2,..., N/2-1, VN>2 (33) 
E|ujv(fc)| 2 - o-^ d/2 , Vife = -N/2,...,N/2- 1, VA>2. (34) 

From this we see that as we refine the grid, limjv^co E|ujv(fc)| 2 = °o for all fc and d > 1. 
Generalizing the discrete version of the iJ s (T d )-norm ([23]) to general dimension d > 1, and 
using the growth rate (p4|) . we find 

JV/2-1 

E||K|||^ = p- d / 2 £ (i + l^l 2 ) s l%(fc)l 2 

fei=-iV/2 
/.jV 

~p- d/4 (l + r 2 Yr d ~ 1 dr~N 2s+d l 2 . (35) 



Jo 

The limits in ([3H|) follow now easily. □ 

In particular, (|3U|) implies that the continuum limit of itjv does not take values in H s (T d ) 
for any s > —d/4, and in particular is not in L 2 (T d ). Furthermore, recalling that the dual 
of H s (T d ) is H- s {T d ), and that C°°(T d ) C H s (T d ) Vs G M, it follows immediately that 

Corollary 4.2. Lei d > 1 and U(ir) = jjir 4 — ^£ 2 - TTiera the stationary solution un of the 
regularized problem \28\) converges in probability to the zero- distribution in the sense of the 
definition given in Coniecture \l.l\ 



4- 2. Simulations 

Similarly to the case of the stochastic heat equation in Section 13.21 we compare 
our results on the decoupled stochastic Allen-Cahn equation now to numerical experiments. 
According to (l34l) the mean of the radial energy density En ([25]) decays like 



EE n (k)~—, k = 1,...,N/2-1, N^oo. (36) 
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Figure 3: 2D decoupled stochastic Allen-Cahn equation . Equation H58H is solved using scheme (I5Q D . 
Parameter values: g = 0.1, cr = 7r/5, T = 2, M = 4 ■ 10 3 . For each N = 2 n (n = 5, . . . , 11), the sample mean 
of Bj\r( K ) over 40 simulations is plotted for k = 1, . . . , N/2. Largest error bar is shown for each N. 

independently of the spectral radius k. This is illustrated by the simulations in Fig.[3l where 
the energy density is plotted for various values of N. Taking into account the loglog-scale of 
the plot, the regular spacing of the Em for increasing N (N = 2 3 ' for j = 5, . . . , 11) illustrates 
the decay in (|36|). A detailed discussion of the employed numerical scheme is presented in 
Section O 

5. 2D Stochastic Allen-Cahn equation 

5.1. Heuristic Considerations 

We turn our attention to full the stochastic Allen-Cahn equation in two dimensions, 
i.e. we consider 



on T 2 , subject to homogenous initial conditions. As outlined at the end of Section [5] it is 
unclear how to make sense of this equation: in fact, we suspect the continuum version of (|37p 
to be ill-posed with respect to any reasonable definition of a weak solution. Therefore, we 
take a different approach and tackle the problem from the perspective of numerical analysis. 
By cutting off high frequency modes in the noise, 



du = [Au + u- u 3 ] dt + adW 



(37) 




(38) 



|fe|<JV 



we obtain the following family of regularized problems 



diiN — [Am at + un — u N ~\ dt + adWu- 



(39) 
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For fixed N, the solution to (|3T))) exists and is unique, but what can we say about the limit 
of un as N — > oo ? In the remainder of this section, we discuss this limit and provide 
evidence for Conjecture 11.11 the main subject of this study. The arguments hereafter will 
invoke elements of the theory of stochastic quantization and we refer the reader to [U for a 
rigorous introduction to the formalism. 

At the end of Section [3] we explained that the term u N in (l39t was expected to be undefined 
in the limit N — > oo. To circumvent this issue, we introduce a constant Cjv > 1, depending 
only on N, and rewrite equation Q39[) as 

<1un = [Aun — {Cn — 1) un — un {u 2 n — Cn)] dt + adWN- (40) 

The key idea is to choose the constant Cn wisely, so that the term un (u n — Cn) approaches 
a well defined limit. More precisely, we choose Cn such that un (u 2 n — Cn) converges to 
: it 3 :, the renormalized cube of u, which we briefly explain here. Given the stationary 
measure fiN of the underlying linear equation 

dv N = [Av N - (Cjv - 1) v N ] dt + adW N , (41) 



x i — y\ x n '. denotes the n-th Wick power with respect to the limit measure limjv-j.oo /Ujv- 
In contrast to limjv_>oo u%, the renormalized cube limjv^oo un {u 2 n — Cn) —■ u 3 : is well- 
behaved in the sense that its generalized Fourier coefficients (: it 3 :,ek) are uniformly L 2 - 
bounded in the probabilistic sense. But what is the correct choice for Cv? In fact, its job 
is to subtract the diverging parts that cause the unmodified cube to be ill-posed, and it is 
given by 

C N = 3E Wv \u n \\ (42) 

where denotes the expected value with respect to jUjv- To determine the measure /ijv, 
we project pXj) onto the Fourier basis to obtain a collection of i.i.d. Ornstein-Uhlenbeck 
processes. From this it easily follows that /xjv is a mean zero Gaussian measure and according 
to (|16l) its covariance operator is given by 

w ,e k if |fc| < N, 



1 if \k\ > N. 1 ' 



Combining (H2l and (jlB")) yields the following equation for Cn 

°n = A E 2 ■ ( 44 ) 

87r2 ,Sv^-l + l fc l 

For all N > 2, equation (j44|) admits a unique positive fixed point, Cn — > oo as N — > oo, 
and the rate of divergence is logarithmic 

C N = ^logjV-t-o(l). (45) 

47T 

To show (|45p , rewrite the sum in f)44[) as a Riemann sum, estimate the latter by an integral 
approximation (for large N), and evaluate the integral to get 



3a 2 




_ 3(j2 ' 




(£) 








47T 



logiV- ^logC/v 



(46) 
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Making now the Ansatz Cn — Ajv log N and plugging it into (jlS) yields the estimate P5j) . 
With this choice of Cn, we go back to the original equation (|40|) and project it onto the 
Fourier basis to obtain 

du N (k) = -(c N -l+ |fc| 2 ) u N (k)dt - {u N (u 2 N - C N ), e k ) dt + a N (k)dp k {t), (47) 

where UAr(fc) :— a if |fc| < N and aN(k) := if |fc| > N. Recalling that Cn — > oo and 
(un(u n — Cjv),efc) — > (: u 3 (which are uniformly L 2 (P)-bounded), the first term in 

(14"T1) dominates the second term as N becomes large. Hence 

du N (k) ~ - (C N - 1 + |fc| 2 ) u N {k)dt + a N (k)d/3 k (t). (48) 

The solution to P5|) is an Ornstein-Uhlenbeck process for |fc| < N and a decaying exponential 
for \k\ > N, and we conclude that for large times 

E\u N (k)\ 2 ~ , ^ (fc) -v, t^oo. (49) 

2 [Cn - 1 + |fc| 2 

From this we deduce two important facts: first, for a fixed cut-off N, the energy in Fourier 
space essentially decays as 

E\u N (k)\ 2 2> A > 1 fixed, \k\ < N. (50) 

|fc| 

Second, for a given mode k, the decay as A — > oo goes as 

E|M*)| 2 ~ tA?' fixedfc - ( 51 ) 
log TV 

Furthermore, given the joint dependence on /c and iV in (|49p . we can estimate the continuum 
limit of the £P(T 2 )-norm of un- For s = 0we get the L 2 (T 2 )-norm 

^ll u Jv||o = ~2 E ^ |u7v(fc)| 2 ~ log f 1 + ~^=^ ~* 00 asA^oo, (52) 
|fc|<iv " V v w/ 

because (jlS"]) implies that N~ 1 ^/Cn — > as N — > oo. On the other hand, if s < 0, we can 
use Lebesgue's dominated convergence theorem to get 

2 

n\uN\\ 2 s = Y E (i + ifci 2 ) s Ei^(fc)i 2 

|fc|<JV 

— (T 2 7r / lrojvi -s t= — dac as JV ^ oo, Vs < 0. (53) 

Jo x + c n 



Note that the results (1521) and ([531) have been derived under the assumption that the field has 
reached the steady-state distribution (t — > oo). However, since the decay rate (Cn — 1 + |fc| ) 
in (|48[) tends to infinity in the continuum limit, the damping occurs infinitely fast as N — > oo, 
and we anticipate the limits (|5"2"|) and (|5"3")) to hold for any finite time t > 0. In other words 
we propose the following conjecture. 



Conjecture 5.1. The solution un to the regularized problem L39\) admits the following 
limit. For all t > 0, 
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An immediate corollary to Conjecture 15. H is Coniecture ll.il un converges in probability to 
the zero-distribution. Even though these two conjectures do not tell us anything about the 
nature of the continuum solution of the stochastic Allen-Calm equation (I57|) . they suggest 
that the approximation sequence (1591) does not admit a meaningful limit: even though the 
limit takes values in the distributional spaces H s (T 2 ) almost surely, for all s < 0, it is a 
trivial limit as it projects every test function onto 0. 

To gain a better understanding of the nature of this pathology, it is instructive to have a 
closer look at the roles played by the three main ingredients of the regularized stochas- 
tic Allen-Cahn equation (1391) : diffusion, nonlinearity, and noise. Considering only the 
noise, the solution to the corresponding equation, du^ = cLWn, converges to the cylin- 
drical Wiener process W as N — > oo. The latter takes values in H S (T 2 ) almost surely, 
for all s < —1. Adding the nonlinearity to the noise, we obtain the decoupled equation 
dux = \un — u n] dt + dWtf of Section 0J In this case, the un converge in probability to a 
distribution u c which is more regular than W: u c € H s (T d ) almost surely, for all s < —1/2; 
however, u c is the zero-distribution. Finally, when adding the Laplacian to obtain the regu- 
larized stochastic Allen-Cahn equation, duN = [Aitjv + un — u%] dt + dWjy, the regularity 
of the limit as N —¥ oo is conjectured to improve: (|54[) implies that un converges in prob- 
ability to a distribution u c with values in H S (T 2 ) almost surely, for all s < 0. At the same 
time, the presence of the Laplacian is not sufficient to render the limit meaningful: u c is 
again the zero-distribution. 

5.2. Simulations 

The results of the previous section are based on heuristic arguments in absence of 
rigorous proofs. To put them on more solid ground, we performed a series of numerical 
experiments. Similarly to the numerical analysis for the simplified equations in Sections 13.21 
and 14.21 we integrated the stochastic Allen-Cahn equation over the time interval [0, 1] and 
computed the radial energy density Em at t = 1 , see Figure [4j A first observation concerns 
the decay of En(k) for fixed N: Figure [4j\ shows that this decay is proportional to ~ 1/k 2 
for large k, which is in agreement with the estimate (|50p. A second observation concerns 
the decay of En(k) for fixed n. In Fig. [4j3 we zoom into the interval k < 16 of Fig. [4j\: 
rather than converging to a finite value, the energy density slowly decays with increasing 
N. In agreement with estimate (|5ip this decay is logarithmically slow, and hence the ex- 
pected convergence to zero is not visible within the range of computationally tractable TV. 
A detailed discussion of the employed numerical scheme is presented in Section 16.31 

A comparison between the three figures - Fig. [5] for the stochastic heat equation, Fig. [3] 
for the decoupled stochastic Allen-Cahn equation and Fig. 2] for the stochastic Allen-Cahn 
equation - suggests that the latter inherits and combines the main features of the two 
simplified versions: the ~ 1/k 2 decay of the Fourier spectrum as dictated by the Laplacian, 
and the decay (to zero) of individual Fourier modes as dictated by the nonlinearity. 

6. Numerics 

To perform the two-dimensional numerical experiments presented in Sections 13. 2[ 
14.21 and 15.21 the periodic domain [— tt, +7r] 2 was discretized into a regular grid of TV 2 points 
(N even), and the interval of integration [0, T] was divided into M subintervals of length 
At = T/M. Depending on the nature of the equation, the time-stepping was performed 
in either direct space for the decoupled Allen-Cahn equation (|27|) . in discrete Fourier space 
for the heat equation (fT3|) . or using both spaces for the Allen-Cahn equation 1|37|) , Note 
that for all presented experiments we had verified convergence in time by direct comparison 
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Figure 4: 2D stochastic Allen-Cahn equation. The stochastic Allcn-Cahn equation l|60|l is solved in 
two space dimensions using scheme 1611 . Parameter values: a = 6.4 ■ 10 -3 , g = 0.5, e = 27r/5, T = 1, 
M = 2 ■ 10 3 . A. For each N = 2 n (n = 5, . . . , 11), the sample mean of Ej^{k) over 40 simulations is plotted 
for k = 1, . . . , N/2. Largest error bar is shown for each N. B. We zoom in on the first 16 modes of panel A. 

For each N = 2™ (n = 5, , 11), the sample mean of Ejv( k ) over 40 simulations is plotted for K = 1, , 16. 

Largest error bar is shown for each N. 



of the solutions obtained with time steps 2At and At, respectively. Before presenting the 
numerical schemes, a comment about diffusion and reaction rates is in order. Note that 
equations (f2Tl) and (1571) are free of constants, except for the noise intensity a. In 

fact, we had eliminated all other constants by rescaling time and space in the corresponding 
evolution equations. For numerical purposes however, it is convenient to use the same 
space and time scales throughout the experiments, and hence we re-introduced diffusion 
and reaction rates (as specified below) in the previously rescaled terms. 

6.1. Stochastic Heat Equation 

Introducing a diffusion constant a > and a decay rate g > 0, the stochastic heat 
equation p^|) becomes 

du = [aAu — gu] dt + adW. (55) 

Linearity of the equation allows us to project it onto the N 2 modes in discrete Fourier space 
as 

du N {k) =-(g + a \k\ 2 )u N {k)dt + adW N (k), (56) 

where Wn is defined in and k — {ki,k%). An implicit-explicit scheme was used for 
the time-stepping, combining the explicit Euler-Maryama scheme [23] for the noise with 
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the implicit trapezoidal scheme for the reaction and diffusion terms. With the notations 
Uj}(k) := iiN(k,t m ) and Ax := 2tt/N, the scheme reads 



uT\k) = u%(k) -^(g + a \k\ 2 ) [u™ +1 (k) + fift(k)] + ^^T, (57) 

where £™ = (77™ + and and {C™} are real i.i.d. Af(Q, 1) random variables. 

This representation of white noise on the grid is dictated by the choice of complex eigen- 
functions {efc} and the fact that the Wiener process W is real. Finally, the (Ax) -1 scaling 
in the noise term of (|57|) can be understood by taking the discrete Fourier transform ([2~Tj) 
of the noise term in discrete direct space, as introduced in (|28|). 



6.2. Decoupled Stochastic Allen- Cahn Equation 

The spatially decoupled stochastic Allen-Cahn equation (|2"T)) is conveniently dis- 
cretized in direct space. Introducing the double-well intensity g > and the notation 
Xj = jj-j, where j = (ji, 32) and ji = —N/2, . . . , N/2 — 1, we obtain the discretized version 

duN(xi) = g [uN(xj) — u 3 y (xj)] dt + adW^ixj). (58) 

Again, an implicit-explicit time stepping scheme was employed for (|58p. combining the 
explicit Euler-Maryama step for nonlinearity and noise with the trapezoidal scheme for the 
linear contribution 



u™+\x 3 ) = u^+g^ [u™ +1 {x 3 ) + u%{x 3 )} - g At{u%(x 3 )] 3 + a^^, (59) 
where {£™} are real i.i.d. AT(0, 1) random variables. 

6.3. Stochastic Allen-Cahn equation 

Since the stochastic Allen-Cahn equation 

du= [aAu + g {u - u 3 )] dt + adW (60) 

involves the characteristic features of both prior models, the diffusion term as well as the 
nonlinearity, there is no natural choice between direct and Fourier space. Instead, we formu- 
late the scheme in Fourier space, go to direct space using the inverse fast Fourier transform 
(IFFT), process the nonlinearity ujy l— ^ u %, and go back to Fourier space using the fast 
Fourier transform (FFT). More precisely, the discretization scheme for (IBTJl) reads 

u™ +l {k) = W) + (9 ~ a \k\ 2 )^ [uT\k) + fijftfc)] - g At ((u^ f , e,)^ + v^-tf, 

X (61) 

where £™ = (77™ -I- with {f?™} and {C™} real i.i.d. Af(0, 1) random variables. In 

(|6T|) we have used the notation 

(x 3 , e k ) N := [FFT [(IFFT [x]) 3 ] ] ^ , (62) 

for x G span {efc : fcx,2 = —N/2, . . . , N/2 — 1}. Straight-forward implementation of the pro- 
cedure (|62l) gives raise to an aliasing error, and we employed the so-called Two-Thirds Rule 
36] to correct for this error. We refer to Appendix A in 17] for a detailed discussion of 
aliasing issues and the Two-Thirds Rule. 
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7. Conclusions and Outlook 

Let us now discuss the implications of our work with respect to the discrepancy of frequent 
use (applied sciences) and virtual neglect (mathematical community) of the white noise- 
driven Allen-Cahn equation in two dimensions. Whereas the well-posedness of the equation 
per se remains an open question, Conjecture 15. II suggests that modeling with the stochastic 
Allen-Cahn equation is indeed problematic: if the mesh size in simulations was shrunk, the 
numerical solutions would converge to the zero-distribution. As illustrated in Fig. [T] there 
is no pattern formation in the continuum limit - the zero-distribution is a trivial solution. 
Therefore, even if one were able to make sense of the continuum equation, the employed 
numerical schemes could not be used to approximate it. In addition to the pseudospectral 
method described in Section 16.31 wc positively tested our conjecture for a finite difference 
method. Since the pathologies do not seem to be caused by subtle deficiencies in the nu- 
merical schemes, but rather by the roughness of white noise itself, we expect Coniecture l5.ll 
to hold for all standard numerical schemes. 

One way to interpret simulations of the stochastic Allen-Cahn equation is to view them as 
numerical approximations of equations driven not by white noise, but by a noise field having 
a finite correlation length A. Not only does this render the model equations well-posed, but 
it can also provide a justification of the numerical experiments: in fact, if A Ax, where 
Ax is the smallest mesh size used in the simulations, the noise is indeed uncorrelated on the 
grid, and the white noise simulations are reasonable. However, such a correction may be 
problematic as A is dictated by the physics of the problem: it might be difficult to determine 
the correct correlation length, or the latter might be too big (A > Ax), in which case the 
white noise simulations are indeed invalid. 

For future endeavors, our work leads to the following questions. First, Conjecture 15.11 poses 
a challenge in terms of conceiving a rigorous proof, and we will address this issue in a 
forthcoming publication [39}. Second, since the irregularity of white noise increases with 
the spatial dimension, we expect equations in three and more dimensions to be ill-posed, 
too. But a generalization of Conjecture 15.11 for d > 3 is yet to be established. Finally, 
the stochastic Allen-Cahn equation is only one specific example of the class of white noise- 
driven nonlinear SPDEs that are used to model physical phenomena in two and three space 
dimensions. It remains to be established whether other models (see Introduction for a 
detailed list) suffer from similar pathologies. 
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9. Appendix 

To prove Theorem 14. II we need the following lemma. 

Lemma 9.1. Let V(x) = jX 4 — ^x 2 be the double-well potential, and let d > 1. Then the 
stationary solution upf of the regularized problem (20) has the following properties: 

(i) EuNixj) = 0, Vji, 2 =-N/2,...,N/2- 1, VA > 2. 
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(ii) E|ujv(^)| 2 ~ aN d / 2 , Vj'i, a = -JV/2,...,JV/2-l, ViV > 2. 

Proof, (i) follows from the symmetry of the potential V. For (ii), we recall that the station- 
ary solution of (1281) has a distribution at every grid point given by 

, , i r 2V(x) i i r (.t 4 /2 - x 2 ) i a^ 2 

2W) = T7 ex P1 f = T7 ex Pi _9 .w/9 r ' ?=7Z2- ( 63 ) 



TV \ ff y/ 2 j A/ \ a 2 // 2 J ' ' 4tt 2 
To simplify calculations, we denote 

— 1 - (27r) V*. (64) 



' • 2o- 2 p d / 2 2ct 2 
Then 

p N {x) = — ■. e c exp {-c{x 2 - l) 2 } 

Ekjv(x-)| 2 = ^ x2pjv ^ rfa: = f™ x2 PN(x)dx = / °°x 2 exp{-c(x 2 - l) 2 } dx = _ 
J r Pn(x)(1x L p N (x)dx J Q exp {— c(x 2 — l) 2 } dx 

(65) 

since the integrands are even. With the change of variables z = \fc-x 2 we get 

R(c) = J™ ^c ex Pi-(z - V~c) 2 } dr* dz 
f °° exp{-(z - Vc) 2 } 2jh7* dz ' 

Now we define 

and use (|64|) and ([66| to see that 

P(c) = y/2ap d '*P(c). (68) 

(2-7r") 2 

Next we establish that P(c) is bounded on c < „ ; -4=> A" > 1. We recall the integral 
representation of the modified Bcssel function K v {z) (47], pl83) 

and use it together with the software Maple to get 



P(6)= 24i)M-^)' (69) 

Exploiting recursive differentiation formulae for Bessel functions |47| . one can deduce from 
(|69|) that P(c) is bounded and decreasing on intervals [0,M], for all M > 0. Recalling (|64|) 
we conclude that P as a function of A" is increasing and uniformly bounded for N > 1. 
Together with JBS]) this yields (ii). □ 
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